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1 Introduction 

The equation of state (EoS) of asymmetric nuclear matter at high densities is a key 
topic for the description of matter inside neutron stars. The determination of the 
properties of asymmetric nuclear matter, such as the symmetry energy {asym) and 
the slope of the symmetry energy (Lq) at saturation density, has been exaustively 
studied in order to better constrain the nuclear matter EoS. However, differently 
from symmetric matter properties that are reasonably constrained, the symmetry 
energy and its slope still large uncertainties in their experimental values. Regarding 
this subject, some studies point towards small values of the slope of the symmetry 
energy [DEIE], while others suggest rather higher values HElEllTllH]. Such a lack 
of agreement raised a certain debate in the scientihe community. In this paper, we 
aim to analyse the role of these properties on the behavior of asymmetric hyperonic 
matter. Using the formalism presented in Ref. [9], which considers many-body forces 
contributions in the meson-baryon coupling, we calculate the EoS of asymmetric 
hyperonic matter (section 2) and apply it to describe hyperonic matter and hyperon 
stars (section 3). 

2 Formalism 

In this work, we use the formalism developed in Ref. [9] that takes into account 
many-body contributions to the nuclear force via a parameterized derivative coupling 
for the mesons. Besides the a and u mesons present in RMF models, we also consider 
the g and 6 mesons, that allow a better description of the properties of asymmetric 
matter, as well as the strange (j) meson, that is important for the description of 
hyperon interactions. Also, as we are foccusing on stiffer EoS, we neglect the strange 
a* meson, since it introduces more attraction to the system. The Lagrangian density 
of the model is (see Ref. [9] for more details): 
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^ - ggbhbQs) - {mb - g^cr - ggbhbSs)] il^b 

+ \ {dixfyd^^a - m^cr^) + ^ {df^dsd^^ds - m|(5|) + ^ + mlto 

+ ^ + mlcl)^(l)^^ + i + mlgl-Qf^ + {id^ - mi) ipi. 

(1) 

The subscripts b and I represent the baryon octet (n, p, A°, S“, S°, S+, S“, H°) 
and lepton (e“, p~) degrees of freedom, respectively. The Dirac Lagrangian density 
for baryons and leptons are in the hrst and last terms, while the other terms represent 
the Lagrangian densities of the mesons: a Klein-Gordon Lagrangian density for the 
scalar a and 6 helds and a Proca Lagrangian density for the vector u, g, and 0 
helds. The baryon isospin component in the 2 ;—direction in represented by hb, which 
correspond to the terms that take the isospin asymmetry into account 

As we use the so called scalar version of the model in this work, the many-body 
forces contributions are introduced only for the scalar meson-baryon coupling con¬ 
stants, present in the effective mass term of Equation [T] The dehnition of the 
meson-baryon couplings are: 

9(7bcr + g5bhbb3\ I 9(TbCr + gsbhbb3\ 

- 7 - 9 ^b, = 1 +--- gsb ( 2 ) 

Qrrib J \ Qmb J 

where the coupling introduces nonlinear contributions, controlled through the C, pa¬ 
rameter of the model and the effective mass of baryons have the following expression: 
= mb- g*ba - glbhbh- 

Also, as we allow hyperonic degrees of freedom to appear at high densities, we 
describe their coupling constants to the uj and 0 vector mesons using SU(6) spin-flavor 
symmetry mm, as follows: 

1 _ 1 _ 1 _ 2^/2 

-guiN — -guiA — -fi'ws — gujE, — -^9ujn — ^g<j>A — 2 ( 7^2 — g<t>Ei (3) 

and also assume an isospin scaling for the coupling to the isovector mesons: 

9eN = l^ggT. = ggE, ggA = 0, gSN = -9513 = fi'5A = 0. (4) 

Finaly, the coupling to the sigma meson is obtained by htting the potential depths 
of the hyperons in nuclear matter [Ill [13] for the values = —28 MeV, = 
-|-30MeV, and = —18 MeV [H]. In what follows, we vary the values of the 

symmetry energy and its slope and analyze their effects on the properties of hyperonic 
matter and hyperon stars. 
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33 
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25 
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25 

90 

13.44 

7.55 
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8.54 

0.059 

0.70 
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33 

100 

13.44 

7.55 

5.84 
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33 
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Table 1; Nuclear matter properties for different parametrizations of the model (dif¬ 
ferent C’s): effective mass of the nucleon, compressibility modulus, symmetry energy, 
symmetry energy slope and (a, a;, g, 6) coupling constants. 


3 Results and Discussion 

Using the mean field approximation, we calculate the EoS of the model from the com¬ 
ponents of the stress-energy tensor, imposing a saturation density of po = 0.15 fm“^, 
and a binding energy per baryon oi B/A = —15.75 MeV. Also, by the fitting of the 
mesons’ coupling constants, we apply the EoS to describe the following nuclear matter 
properties at saturation: effective mass of the nucleon compressibility modulus 
Kq, symmetry energy and symmetry energy slope Lq. These quantities and 

the respective values of the coupling constants are displayed in Table [H for different 
parametrizations. 

In addition, there is still debate in the literature regarding the values of the asym¬ 
metric matter properties at saturation na. For this reason, we vary their values in 
ranges of between 25 — 33 MeV [5] and Lq between 73 — 100 MeV [H UHl HTJ HE] , 
in order to verify their effects on the behavior of hyperonic matter at high densities. 

As discussed in the previous section, the effective mass of baryons depends on 
the ^-parameter which, consequently, is also affected by the nonlinear contributions 
that come from the many-body forces coupling. Hence, each parametrization has 
different particle populations, as one can see in Ref. [9]. Furthermore, as the isospin 
asymmetric properties and Lq) have direct impact on the coupling constants of 
the isovector mesons (p, 5), they also have effect on the particles population. 

It is important to emphasize that, for a hxed value of the symmetry energy, the 
lowest values of symmetry energy slope depends on the many-body forces parameter: 
the smaller the ( parameter, the larger the symmetry energy slope [9]. In particular, 
the smallest values of the symmetry energy slope for the parametrizations used in 
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Figure 1: Particle population dependence on the symmetry energy 0 °^^, for the 
parametrization ( = 0.040 and fixed symmetry energy slope Lq = 100 MeV. The left 
panel shows the population for = 27 MeV and the right for = 33 MeV. 
The x-axis represents the baryon density and the y-axis the fraction of each particle 
species normalized by the total baryon density. 


this work are those presented in Table [T], for each value of symmetry energy. 

Fig. [Hand Fig.[2]show the particle population dependence on the symmetry energy 
and its slope, respectively. In Fig. [H we show the results for the parametrizations 
C = 0.040 and Lq = 100 MeV, for a^y^ = 27 MeV (left panel) and a^ym = 33 MeV 
(right panel). In Fig. [2l we show the results for the parametrizations ( = 0.054 and 
a^ym = 25 MeV, for Lq = 73 MeV (left panel) and Lq = 90 MeV (right panel). 
The results show that lower values of the symmetry energy and higher values of 
the symmetry energy slope favor (disfavor) the appearence of particles of negative 
(positive) isospin. Both quantities {a^y^ and Lq) have proportional effect on the 
coupling constants of the g and 6 mesons, i.e. lower values of a^ym and Lq provide 
weaker coupling of the baryon to these mesons. These coupling constants, together 
with the isospin projection of the particles, play a role on the chemical potential and 
effective masses of baryons calculations and, consequently, on the particle populations. 

Regarding the fraction of strangeness fs, the left panel on Fig. [3] shows that both 
the many-body forces (^-parameter, as well as the slope of the symmetry energy have 
impact on hyperonic matter. On the right panel of Fig. |3l one can also check that 
the effect of the symmetry energy on the fraction of strangeness is negligible, for the 
range of values studied in this work. 

Finaly, we apply our EoS to describe hyperon stars, solving the Tolman-Oppenheimer- 
Volkoff equations, and perform the same analysis as above for the Mass-radius dia¬ 
gram in Fig. m On the left panel of Fig. 01 we show the results for three parametriza- 
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Figure 2: Particle population dependence on the symmetry energy slope Lq) for the 
parametrization ( = 0.054 and hxed symmetry energy = 25 MeV. The left panel 
shows the population for Lq = 73 MeV and the right for Lq = 90 MeV. The axes are 
the same as Fig. [1] 
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Figure 3: Fraction of strangeness fs as a function of baryonic density. The left 
panel shows the results for three different parametrizations = 0.040, 0.049, 0.059) 
and an interval of roughly 15 MeV for the symmetry energy slope Lq, and hxed 
symmetry energy = 29 MeV. The right panel shows the strangeness results for 
the parametrization ( = 0.040 and symmetry energy slope Lq = 100 MeV. 
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Figure 4: Mass-radius relation for different values of symmetry energy and its 
slope. The left panel shows the results for three different parametrizations = 
0.040, 0.049, 0.059) and an interval of roughly 15 MeV for the symmetry energy slope 
Lq, and fixed symmetry energy = 25 MeV. The right panel shows the strangeness 
results for the parametrization ( = 0.040 and symmetry energy slope Lq = 100 MeV. 


tions (C = 0.040, 0.049, 0.059) and different values of Lq, for = 29 MeV. The 
results for the impact of the symmetry energy on the Mass-radius diagram are shown 
on the right panel of Fig. 01 for the parametrization ( = 0.040 and Lq = 100. The 
results show that the symmetry energy and its slope have small effects on the max¬ 
imum masses of the stars but, on the other hand, have signihcant impact on the 
radii. In particular, a shift of 15 MeV on the slope of the symmetry energy Lq 
changes the radii of the canonical star (1.4 Mq star) in roughly 0.5 km. Similarly, 
a shift of 6 MeV on the symmetry energy a^y^, changes the radii of the canonical 
star in 0.4 km. Also, the many body forces parameter presents strong influence on 
both the maximum mass and radius of the stars, with the highest maximum mass 
{Mmax = 2 .I 5 M 0 ) being provided by the parametrization ( = 0.040 [9]. 

In this work we have analyzed the effects of many-body forces, symmetry energy 
and its slope on the properties of hyperonic matter and hyperon stars. We have 
verified that the many-body forces contributions play an important role on the particle 
populations and on the macroscopic properties of hyperon stars. Also, we showed 
that the asymmetric matter properties have non-negligigle influence on the particle 
populations and strangeness distribution of such systems, as well as have impact on 
the determination of the radii of hyperon stars, as already suggested by several studies 

[MUnillH]. 
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